NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



NASA Technical Memorandum 81201 

(NASA-TM-81201) MODIFIED ITERATIVE EXTENDED 
HUCKEL. 2: APPLICATION TO THE INTERACTION 

OF Na (♦) , Na(+)(ag.), Mg (+) -2 (ag. ) WITH 
ADENINE AND THYMINE (NASA) 28 p 

HC A03/MF A01 CSCL 2QH G3/72 


1/8 0-25109 


Unclas 

17902 


Modified Iterative Extended Huckel: 

II. Application to the Interaction of 
Na + , Na + (aq.), Mg + 2(aq.) 

With Adenine and Thyamine 

S. Aronowitz, R. MacElroy, and S. Chang ; 

r 


April 1980 


NASA 

National Aeronautics and 
Space Administration 



NASA Technical Memorandum 81201 


Modified Iterative Extended Hiickel: 

II. Application to the Interaction of 
Na + , Na+(aq.), Mg+2(aq.) 

With Adenine and Thyamine 

S. Aronowitz 

R. MacElroy 

S. Chang, Ames Research Center, Moffett Field, California 


l\J/\SA 

National Aeronautics and 
Space Administration 

Ames Research Center 

Moffett hold California 0403ft 


Modified Iterative Extended HUckel: II. Application to the 

i i i n t j [ /^ 

Interaction of Na , Na (aq.), Mg , Mg (aq.) with Adenine and Thymine 

S. Aronowitz,* R, MacElroy, and S. Chang 

Extraterrestrial Research Division, Ames Research Center, 

NASA, Moffett Field, California 94035, U.S.A. 

ABSTRACT 

Modified Iterative Extended HUckel, which includes explicit effective 
internuclear and electronic Interactions, is applied to the study of the 
energetics of Na , Mg , Na (aqueous), and Mg (aqueous) Ions approach- 
ing various possible binding sites on adenine and thymine. Results for 
the adenine + ion and thymine + ion are in good qualitative agreement 
with a limited basis set ab initio work on analogous systems (adenine 
and uracil ) 5 energy differences between competing sites are in excellent 
agreement. Hydration appears to be a critical factor in determining 
favorable binding sites. The result that the adenine N1 and N3 sites 
cannot displace a water molecule from the hydrated cation (Ion*4H 2 0) 
indicates that they are not favorable binding sites in aqueous media. 

Of those sites investigated, 04 was the most favorable binding site on 

4 - 

the thymine for the bare Na (a full determination was not made with 
4*2 

Mg ). However, the 02 site was the most favorable binding site for 

either hydrated cation. Adenine and thymine results are consistent with 

+2 

experimental findings [Mg (aqueous)]. 

National Research Council Associate. 
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1 * In^dujptlo n 

An ah e inltJLc> investigation of the 1 Interaction between cations and mole- 
cules of biological interest has been pursued In a series of recent 
articles [1,2,3*43 • The motivation, in part, arises from the fact that 
standard semiempirlcal methods either fall completely or are inadequate 
when applied to systems which are sensitive to charge transfer. For 
example, the failure of C.ND0 to correctly predict the most likely bind- 
ing site of a cation on the uracil molecule [2] Is due to its inability 
to deal properly with charge transfer. Similarly, Iterative Extended 
HUekel (hereafter, IEH) overestimates charge transfer and exhibits prob- 
lems in achieving self-consistent charge distributions. These fallings 
are discussed in paper I [5] of this series. Modifications to I Ell, 
which include explicit effective internucloar and electronic interaction 
terms, are detailed in paper I. These modifications should compensate 
for some of the inherent deficiencies of IEH, such as the problems 
associated with charge transfer. Therefore, the application of the 
modified IEH program to catlon-hiomolecular systems is a critical test 
of the usefulness of the modifications to IEH. Adenine and thymine were 
chosen for study as substrates for ion-binding because they are rela- 
tively simple; moreover, calculations have been carried out either with 
these or similar structures [2,4]. The adenine and thymine struc- 
tures [6], along with the possible binding sites that were considered, 
are depicted in Figs-. 1 and 2. The possible binding sites were limited 
to those considered in Refs. [2,4]. Thus, N9 on adenine, or Ml, N3 on 
thymine were not examined* The workers in Refs. [2,4] did not attempt 
a global search; coirs equently, since our interest was to make as direct 


- 3 ~ 


a comparison as was feasible, we also restricted our probing of the 
potential surface to outs within the plane of the adenine and thymine 

4" *j*o 

molecules. The hare Na and Mg “ cations wore used to probe those sites 
noted in Figs. 1 and 2, After determining the most favorable sites for 
the hare ions. Imitated forms were then introduced and the most favor- 
able sites redetermined. This more complete modelling of catlon- 
Momoleeule binding lias not been attempted previously, to our knowledge. 
The modified JF.H program has been described in detail in paper I; how- 
ever, a brief description will be offered in the following section. 


2 * Modified lUtrat l ve Extended i Hickel 


The program builds upon and extends convent ional TE11. Explicit effec- 
tive nuclear interaction terms are included. Let 7,^ (valence) be the 
resultant nuclear charge after assuming that the core electrons com- 
pletely shield their respective number of positive charges. Between 
the ith and jt.h atoms, which are separated a distance R, we set 


(?. 7, ) 

{ i ,] effective 


ZjCvnloncc) - £ l’ k(1) W.-Wp^, 
icd) 

Zj (val mC0 ) - £ 
m(j) 


( 2 . 1 ) 


Here, (0,-aR) is the probability of finding an electron in the kth 

orbital located on the lfch atom in a spherical volume of radius ^R. 
Similarly defined is ^ (O^R) . The term is the electronic 

density assigned to the kth orbital after the JLth iteration. The term 

U) 


P 


m(j) 


lias a similar significance with regard to the mtli orbital on the 


jth atom. The effective product 
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(Z iV effective * Z i( vnlence > z .j (valence) , R + 0; 


(2. 2d) 


on the other hand, 


effective * 




(2.2b) 


where q. 0 ,q. 0 are the net charges on the ith and jth atomic centers. 
The contribution to the total energy made by the effective nuclear- 
nuclear interactions after the ith iteration is given by 


E w 

effective nuc-nuc 


S (Z iV effective /R ij 

i^j 


( 2 . 3 ) 


in atomic units. In a similar vein, the molecular environment is sensed 
more directly by introducing coulombic modifications to the valence 
orbital ionization energies. Here, an electron in the kth orbital on 
the ith atom (which also is taken to be the ith atomic type, for con- 
venience) interacts with the electronic densities assigned to the 
orbitals on the X centers. A set of variables, which depends only on 
atomic types, is formed, and the variables are determined either by 
requiring specific charges (as was done in paper I) or are determined 
by variational considerations. This increment to the valence orbital 
ionization energy is expressed as 


k(i) 


j k(l),3 


fid 


Ep 


U) 

s(j) 






The (a..) are the set of variables depending only on atomic type; 
p^' ^ is the electronic density assigned to the s. orbital on the j_th 
atom after the &th Iteration. The coulombic integral Involves the two 


atomic centers, j. and j., through their respective orbitals, and 

^s(j)' ^ ie * n * t * a -** diagonal elements of the Hamiltonian matrix are now 
given by 


1! k(l),k(l) “ ll k(i),k(i) (0) + h k(i) ’ 


(2.5a) 


The off-diagonal elements 


H 


k(i),s(j) * ,sKS k(i),s(j)[ H k<i), 


k(i) + H s(j) 


»s(j)] ’ 


(2.5b) 


The symbol K is an interaction constant set equal to 1.89; 
is the Integrated overlap between orbitals k on the 1th atom, and £ 
on the atom. The eigenvalues of the Hamiltonian matrix may now con- 
tain two-electron contributions, but these are Included in a fashion 
which permits the energies to be treated as though they were one-electron 
energies. The total molecular energy after the ith Iteration for a sys- 
tem composed of m atoms and N orbitals is given by 


E 


U) 

total 


N m 

Lvi l) + EE( 2 i Z J )e£L C iv a /li « < 2 - 6) 

u— 1 i^j 


(0 ) 

where n u is the occupational number and * is the energy of the utli 
molecular orbital. The first term on the right-hand side of Eq. (2.6) 
obeys a variational principle, from paper I. The E tota ]_ does not, in 
general, since some of the (Z.^ ^ effective terms ma y $ row more posi- 
tive with successive iterations, before charge self-consistency is 
achieved. Standard XEH parameters were used in the calculations, and 
these are collected in Table 1, No attempt was made to optimize these 


parameters for the systems considered. Within the framework of the 
Modified IEH, however, the energy values do represent optimum values 
given the constraints of stability (oscillation) and minimum electronic 


energies. The coefficients a^ in Eq. (2.4) were determined to ±0.005. 
All calculations were performed on the NASA-Ames CDC 7600. 


3. Results 

3. 1 Adenine 

+ +2 

Adenine + ion results for both Na and Mg are displayed irt Fig. 3. 

These ions were restricted to pathways in the plane of the molecule which 

bisected the bond angle of a site. Energy values are relative to the 

isolated adenine, i.e., dE =■ E (adeniiie + lon) -E adenlnB , since the valence 

shell for the Isolated Na , Mg ions are empty [E^ Qn (valence electrons) 

- 0], Net charge on the sodium never exceeds +1. All calculated points 

4 - 

for the adenine + Na thus fall within the range of the Hamiltonian 

matrix algorithm. The algorithm adopts a linear extrapolation procedure 

to determine the valence orbital ionization energies when the atomic net 

charges range from -1 to +1. The algorithm is no longer reliable when 

the magnitude of the net charge |q| >1. We chose not to expand the 

range of reliability. Because of this, the bare Mg ion case presented 

some difficulties. The net charge for Mg along the solid line portion 

in Fig. 3 was < +1. Over the dashed portion of the Mg-Nl curve, the net 

charge is greater than 1.1, and at 4.0 A it equals +1.8. The energies 

of the Mg-Nl curve at 2.5 A, 3.0 A, and 3.5 A were fitted to a quadratic 

2 

equation form, y = Ax + Bx + C, in order to determine the location of 
the minimum. Cramer's rule was used to solve the resulting set of 
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simultaneous equations. The extrapolated minimum and energy are 3.93 A 
and -25.87 eV, respectively. The energy is smaller in magnitude than the 
energy associated with the Mg-N3 minimum. Therefore, the N3 site is pre- 
dicted to be the favorable one for Mg +2 . 

The preferred interaction site for the bare Na is Nl, which is in 

accord with Ref. [4]. The Mg case was not considered by these workers, 

* 4 * 

The difference in energies between Na-Nl and Na-N3 when the Na is at 
either minimum is roughly 0.08 eV (1.8 kcal/mole). This energy differ- 
ence is consistent with the behavior of the Na-Nl and Na-N3 curves over 
the range of 2.5 A to 3.5 A. Reference [4] displays a difference of 
2.4 kcal/mole. Our equilibrium distance, about 3.5 A, is greater than 

the 2.15 A value adopted in [4]. Energy differences between sites and 

+ 

general patterns for the adenine + Na system agree remarkably well 
with [4]. 

The hydrated ion complex was constructed in the following manner: Four 

water molecules were tetrahedrally disposed about the ion (Na or Mg) 
with the oxygens placed at the vertices of the tetrahedron. The tetra- 
hedral structure was chosen because it gives a "saturated" geometrical 
shape which is consistent with the experimentally deduced solvation 
numbers for Na and Mg [7]. Calculations involved only the simultaneous 
and equal dilation or contraction of distances between the ion and the 
four water molecules; that is, the tetrahedron was uniformly expanded 
or reduced in size. The equilibrium distance between the oxygen of the 

water molecules and the ion (at the center of the tetrahedron) was found 
0 

to equal 2.4 A for both Mg and Na using the modified XER. This distance 


Is In good agreement with fits for Mg [7] and Is In reasonable agreement 

4 * 

with the ab inltlo calculations for Na and a single water molecule £8 ] . 

It is to be expected that some expansion in the solvation shell will 
occur when Na* Is surrounded by four waters. One water was removed for 
the Interaction between a hydrated ion and the adenine or thymine. The 
oxygens of the remaining three waters formed an equilateral triangle 
which served as the base of a triangular pyramid whose vertex was at the 
adenine N1 or N3, for example. Figures 4 and 5 typify the approach. 

The ion is in the plane of the adenine. The pathways bisect the angles 
(C2-N1-C6) and (C4-N3-C2). The water hydrogens were disposed so as to 
reduce Interactions between various pairs, while also minimizing inter- 
action between these pairs and the adenine or thymine. The adenine 

resultB are displayed In Fig. 6. The net charge on the Mg was less than 

•4*2 4* 

+1 at all distances considered. In both cases, (Mg* 31^0) and (Na * 311,0) , 

the N1 site is -avored. In fact, the Mg'SHgO-NS curve Is repulsive. The 

binding energies, 4E blnd - E adenlne , nre 

considerably reduced with respect to the adenlne*ion binding energies. 

The pathway examined for Mg'SHgO-Nl clearly is not an optimal one when 
o 4*2 

R > 3.9 A, since the isolated (Mg'SH^O) and isolated adenine system is 
more stable than the interactive one. The most important result Is 
obtained when the displacement energies are combined with the binding 
energies. That in, we ask whether the adenine N1 or N3 nitrogens can 
displace a water from the Ion* 41^0 complex. The total energy difference 
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^displacement + binding " E ion*3H 2 Q + ^0 " E ion*4H 2 0 + AE bind* 

“ E H 2 0 ” E ion*4H 2 0 + doenine* Ion *31120 
ddenine * 


The results are compiled In Table 2 and show that AE Is positive for 
either Mg or Na, Treated as an equilibrium problem, the overall dis- 
placement plus binding reaction for Nl, for example, has an "apparent M 
equilibrium constant 

(adenine • ion* 3H 2 0 ) n ^(H 2 0) 

K(Nl) “ (adenine) (lon*4H 2 0) ~ ex ^ ^displacement -t binding/^) ’ 

(3.2) 

At room temperature (300 K), the approximate value for this equilibrium 
constant when the reaction Involves the hydrated Na Ion equals 

K(N1) * 4,0 x IQ" 31 , 

while the equilibrium constant, when the hydrated Mg Ion Is involved, 
equals 

K(N1) e 2.5 X 10“ 45 . 

Similar results are obtained with the hydrated Na ion for K(N3) ; however, 

4*2 

since the magnesium complex (Mg*3H 2 0) does not bind with the N3 site, 
the equilibrium constant K(N3) = 0, for all temperatures. Therefore, 

Nl or N3 will not displace a water molecule from the complex lon*4H 2 0; 
i.e., the displacement plus binding reaction does not oc r to any 
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x a.0 

significant degree. Hence, we expect that neither Na nor Mg will 
bind to the adenine N1 or N3 in an aqueous medium. There appears to be 
experimental confirmation with regard to the Mg (aqueous) cose [9]. 

3, 2 Thymine 

JL 

The results for thymine + ion are displayed in Fig. 7. Only the Na 

results for approaches to 02 or 04 are presented. The Mg case is not 

displayed because the net charge on the magnesium was greater than +1.5 

at all distances for either an 02 or 04 pathway. The energy algorithm, 

as noted in section 3A, breaks down when the net charge on an atom 

+ 

exceeds +1. In the Na case, the 04 site is more stable than the 02 
site by ca. 10 kcal/mole. Reference [4] finds the 04 site in uracil 
more stable than the 02 site by 4.2 kcal/mole. Binding to the ethylenic 
double bond (C5-C6) is highly unfavorable from our calculations. The 
curves labelled Mg-C ■ C, Na-C » C are repulsive. This is in agreement 
with Ref. [2] with respect to Mg . Energy differences (e.g., 
%a-04’"^'Na-02^ a S cee more closely with Ref. [2] than do the separate 
magnitudes, which mirror the behavior found in the adenine case. The 
results are reversed when the ions are hydrated (see Fig. 8). That is, 
the 02 site becomes the favored site. This arises because of interac- 
tions between the waters and the methyl group attached to C5. Note that 
when the Na*3H 2 0 and Mg* 3^0 entities are sufficiently distant from 
the 02 and 04 sites (>3.0 & for Na*3H 2 0 and >3.4 A for Mg*3H 2 0), the 04 
pathway is favored. It is also noteworthy that neither 02 nor 04 path- 
ways are favorable for Na*3H 2 0 at distances greater than 3.4 A. Assuming 
that these repulsive regions are not artifacts of the method, it then 
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follows that It may not suffice to determine binding energies when eval- 
uating the ease with which binding can bo expected to occur, since the 

repulsive regions favor the recombination of (Na* 3H 2 Q) + with H 2 0 to form 
4 * 

(Na-AHgO) . As with adenine, displacement and binding energies were 
calculated. Results are presented In Table 2. Both sites, 02 and 04, 
exhibit + dlsplocemen| .-0. Thus, the th^ine can displace n 

water from the lon*4H 2 0 complex. The prediction that 02 Is a favored 
site appears consistent with experimental evidence for Mg(nqueous) [9]. 

4 . Discussion and Conclusion 

4* 4* 4-2 

Our results for adenine + Na and thymine + (Na , Mg ) are In good 
qualitative agreement with ab Initio model work on adenine [4] and 
uracil [2]. Energy differences between competing sites, where compari- 
sons can be made, are In excellent agreement, Absolute energy magni- 
tudes are not in agreement; these magnitudes roughly correspond to what 
one might expect from CNDO. However, CNDO does fail when applied to 
uracil (see Ref. [2]), while the Modified IEH results are In accord with 
the conclusions drawn from the ab lnltlo work. Thymine does differ from 
uracil In that a methyl group replaces the C5 hydrogen; no differences, 
with regard to sterlc effects, are to be expected with the bare Ion 
interacting with either the thymine or uracil. The water molecules of 
the hydrated Ions, on the other hand, do interact with the methyl group; 
these interactions lead to 02 being the most favorable site rather than 
04. 

No attempt was made to optimize any of the standard IEH parameters in 
the systems studied. It is indeed gratifying that the results agree to 
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the extent that they do with both ab Initio work on bare Iona and 
experimental trends for hydrated lone. The Modified IEH allows inex- 
pensive calculations to be performed on systems where standard IEH can- 
not be used at all, though experience with standard IEH remains a 
valuable guideline. Thus, we feel that the results obtained in this 
communication are sufficiently encouraging that one could hope to model 
other systems of experimental, Interest which might provide bridges to 
span the gap between theory and experiment. 
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our attention and Dr. Jack Egan for the graphics work (Figs. 4 and 5) 
using the Ames Interactive Modeling System. 
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Table 1. 

Standard Extended ! 

UUckel parameters for 11, G, N, 0, 

No, Mg 


A. Basis 

sot exponentials 

(Bohr radii) 



11 

5 Is 

- 1.0000 




C n 

C 2s 

- 1.6083 

5 2p - 1.5679 



N U 

r, 2s 

- 1.9237 

C 2 p ■ 1.9170 



n a 

r, 2s 

* 2.2458 

r, 2 - 2.2266 



Na a 

C 3s 

■ 0.8358 

r, 3p * °’ 8358 



Mg® 

r> 3s 

- 1,1025 

? 3 - 1.1025 


B, Valence orbital ionisation energies, eV 


M + ■> Mi + (e) 


M -> M + + (e) 

M + (e) M“ 


ii b 

-15.36 


-13.60 

-0.76 

(Is) 

c b 

-31.65 


-20.78 

-8.89 

(2s) 


-23.61 


-11.32 

-0.87 

(2p) 

N b 

-39.99 


-25.16 

-13.80 

(2s) 


-29.65 


-14,16 

-1.08 

(2p) 

0 C 

-47.84 


-32.37 

-19 . 24 

(2s) 


-33.63 


-15,86 

-1.93 

(2p) 

m 1 ' 

-15.15 


-5.09 

+0,96 

(3s) 


-10.74 


-2.96 

+1.81 

(3p) 

Mg d 

-18.92 


-7.60 

+0.47 

(3s) 


-14.29 


-4.47 

+1.40 

(3p) 


a Re£. [10], Table I. 
b Re£. [11], Table. 5. 
c Ref . [12], Table I. 

^Adopted from Ref. [13], Tables D.l, D.4. 
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Table 2. Site binding and displacement energies for adenine and thymine 

+ +2 

interacting with hydrated Na » Mg ions 


Complex 

Site 

^equilibrium ^ 
:Lon-slte ' 

AE1> eV 
(binding) 

AE2^ eV 
(binding + 
displacement) 

Adenine* Na • 3H 2 0 

N1 

3.0 C 

-0.78 

+1.82 


N3 

3. 0 C 

-0.65 

+1.96 

Adenine*Mg • 3H 2 0 

N1 

3.5 C 

-0.26 

+2.67 


N3 

- — 

NO BINDING 

— 

Thymine *Na*3H 2 0 

02 

2. 5 C 

-3.26 

-0.65 


04 

3. 0 C 

-2.81 

-0.21 

Thymlne*Mg* 3H 2 0 

02 

2.5° 

-6.11 

-3.17 


04 

2. 5 C 

-4.67 

-1.74 


a AEl = E(adenine, thymine* ion *3H 2 0) - E(adenine, thymine) - E(ion*3H 2 0). 
b AE2 = A El + E(ion*3H 2 0) + E(H 2 0) - E(ion*4H 2 0). 

c These values are reliable (within the context of the calculations) only 
to ± 0.3 A. 
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Fig. 1 

fir. 2 

Fig. 3 
Fig. 4 

Fig. 5 
Fig. 6 


FIGURE TITLES 


. The adenine molecule 

The thymine molecule 

. Potential curves (AE vs. R) for the bare Na , Mg interacting 
with the Nl, N3, N7 sites on adenine. The cations are brought 
in along paths which are in the plane of the adenine and which 
bisect the respective site angles C2-N1--C6, G2-N3-C4, C5-N7-G8, 
See text for description of dashed line 

. Hydrated cation (three water molecules) approaching the Nl site 

on adenine. The oxygens of the water molecules are. at the 

vertices of the equilateral triangle which forms the base, of 

a pyramid Involving Nl. The apex of the pyramid is occupied 

by the site. One water is above the plane of the paper. The 

0 

other two waters are below. The cation is 4.0 A from the Nl 

. Hydrated cation (three water molecules) approaching the Nl 
site on adenine. Deployment of the waters is the same as in 
Fig. 4 except that the distance between the cation and the Nl 
site is 2.5 A 

. Potential curves (AE vs. R) for hydrated Na , Mg approaching 
the Nl, N3 sites on adenine. The paths are the same as those 


in the bare ion case 
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FIGURE TITLES, continued 

Fig. 7. Potential curves (AE vs. R) for the bare Na , Mg approaching 
the 02, 04 and the ethylenie double bond (C5-C6 in Fig. 2) 
sites on thymine. Paths are in the plane of the thymine and 
are at 180 deg. with respect to C2-02, C4-04. The path to the 
ethylenie double bond disects the bond. The numbers in paren- 
theses along the Na-04 curve is the calculated net charges on 
the Na 

Fig. 8. Potential curves (AE vs. R) for the hydrated (three water 

4* 4*2 

molecules) Na , Mg approaching the 02 and 04 sites on thymine. 
Paths are identical to the bare ion case 
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